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Abstract. We address the problem of computing temperature correlation func- 
tions of the XXZ chain, within the approach developed in our previous works. 
In this paper we calculate the expected values of a fermionic basis of quasi-local 
operators, in the infinite volume limit while keeping the Matsubara (or Trotter) 
direction finite. The result is expressed in terms of two basic quantities: a ratio 
p{C) of transfer matrix eigenvalues, and a nearest neighbour correlator ^). We 
explain that the latter is interpreted as the canonical second kind differential in 
the theory of deformed Abelian integrals. 



1. Introduction 

The present article is a continuation of the paper [2j , which was written almost a 
year ago and was dedicated to the memory of Alyosha Zamolodchikov. It so happens 
that the topic we discuss this time is not too far from a domain in which he made 
giant footsteps. So, life goes on, but there stays a painful sorrow caused by his early 
death. 

Consider the XXZ spin chain with the Hamiltonian 

oo 

(1.1) H=l Z H^Ui + 44^1 + , A = i(g + g-i), 

A;=— oo 

where a" (a = 1,2,3) are the Pauli matrices. To avoid technicalities, in this Intro- 

oo 

duction let us accept fll.ip as a formal object acting on = 'S) We shall 

j=~oo 

touch upon the limit from a finite chain in the body of the text. In the papers [1], 
[2], we studied the vacuum expectation values (VEVs) 

(vac|g2.5(0)|^ac) " 

Here |vac) denotes the ground state eigenvector, S{k) = | Yl'j=-oo '^j' is a local 
operator. We have obtained a description of (11.21) in terms of fermionic operators. 
For that purpose, it was essential to consider operators of the form q'^^'^^'^^O, which 
we call quasi- local operators with tail a. 

An important generalisation of our results was proposed by Boos, Gohmann, 
Kliimper and Suzuki [1]. They gave evidences that our fermionic description works 
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equally well in the presence of a finite temperature and a non-zero magnetic field: 



;i.3) (g2"5(0)Q^ 



XXZ, f3,h 



Trs (^e'^^+'^Sq2aSio)j 



where Trg stands for the trace on ^s- For f3 oo and h = 0, the expectation 
value (11.31) reduces to (11. 2p . For us this was quite an exciting development, because 
it shows that the fermionic structure is not a peculiarity of VEVs, but is rather 
a reflection of a symmetry hidden deep in the model. It should be said that in 
the paper [1] the expectation values (II. Sp were not considered in full generality. 
The formula expressing them in terms of fermionic operators was formulated as a 
conjecture, which was checked in some particular cases but was left unproved. 

The first question which we asked ourselves was, why not to add other local inte- 
grals of motion to —(3H + hS in (11.30 . The physical meaning of such a generalisation 
is obscure, but it should be possible for integrable models. This question, together 
with an intuition coming from the papers [6], [9], led to the following generalisation 
of ( 11. 3p . Along with the space i^s, consider the Matsubara space i^M, 

(1.4) = C^'^"^' ® ■ ■ ■ ® C2'"+\ 

with an arbitrary spin Sm and a spectral parameter Tm attached to each component. 
The generalisation of (11.30 is given by the following linear functional 



, ^ . TrsTrM(Ts,Mg^'^^+^"^WO 
1.5) Z'^|g2"^Wo}- ^ 



TrsTrM(Ts,Mg2-s+2«s(0)j 

Here Ts,m. denotes the monodromy matrix associated with i^s ® (see ( 12. 2p ). 

The idea behind the generalisation ( 11.50 is simple: for whichever spins and that 
we put in the Matusbara direction, TrM(rs,M) commutes with Hxxz- One expects 
that using cleverly this arbitrariness in the definition of ^m, it should be possible to 
reproduce any function of local integrals of motion under the trace. In particular, 
in order to reproduce (11.31) from (II. 5p . one has to take special inhomogeneities and 
then to consider the limit n —>■ oo. This point is explained in detail in [6], [9]. In 
the present paper we compute Z'^ for finite n, leaving the discussion of the limit 
for future publication. We would like to emphasise, however, that this limit is not 
complicated. For finite n, Z'^ will be expressed in terms of only two functions, p{C), 
a;(C,0 (I1.12P below) and one needs only to take the limit of them. Let us 
explain all that in some more details, starting from our fermionic operators. 

For the moment we forget about the Matsubara direction, and concentrate on the 
description of the operators acting on i^s- The logic of our papers PP, [2] is close to 
that of CFT: we describe the space of quasi-local operators as a module created from 
the primary field g^^-^C^) by creation operators. We recall below the main features 
of the construction in [2]. 

We say that X = g2aS'(o)Q jg quasi-local operator with tail a if it stabilises 

3 

outside some finite interval of the infinite chain: to q°"^J on the left and to Ij on the 
right. The minimal interval with this property is called the support of X. The spin 
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of X is the eigenvalue of S(-) = [5, ■] where S = S{oo) is the total spin operator. 
We denote by W„ the space of quasi-local operators with tail a, and by its 
subspace of operators of spin s G Z. Consider the space 



a—s,s 

— oo 



On this space we defined the creation operators t*((^), b*((^), c*(C) and annihilation 
operators h{(), c{(). These are one-parameter families of operators of the form 

oo 

p=i 

oo oo 

p=i p=i 

oo oo 

b(c) = - i)"'b, , c(c) = c - 1)"% • 

p=0 p=0 

The operator t*((^) is in the center of our algebra of creation-annihilation operators, 

[t*(Cl),t*(C2)] = [t*(Cl),C*(C2)] = [t*(Cl),b*(C2)] = 0, 

[t*(Ci),c(C2)] = [t*(Ci),b(C2)] = o. 

The rest of the operators b, c, b*, c* are fermionic. The only non- vanishing anti- 
commutators are 

[b(Cl), b*(C2)]+ = -^(C2/Cl, «) , [c(Cl), C*(C2)]+ = ^(Cl/C2, a) , 

where 

Each Fourier mode has the block structure 



1.7) t* : W,_,,, ^ W, 



a—s,s 



h*, Cp : Wa-s+l,s-l y^a-s,s , C*, bp : Wa_s-l,s+l y^a-s,s ■ 

Among them, r = plays a special role. It is the right shift by one site along 
the chain. Consider the set of operators 

(1-8) ^'"t;---t;,b;---b;c;...c;(g2"^(°)), 

where m e Z, j, G Z>o, Pi > ■ ■ ■ > Pj > 2, gi > • ■ ■ > > 1 and ri > • • ■ > > 1. 
It can be shown that (11. 8p constitutes a basis of V^a,o (we postpone the proof to 
other publication). 
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Now we start to consider the spaces i^s and together. We shall prove that 

(1.9) Z'^{t*(C)(X)}=2p(C)Zn^}, 

(1.10) z^{h*{0{x)} = ^/^(C,0^ic(0(x)}^, 

r 

(1.11) z-{c*{c){x)} = /^(e, Oz^{HO{x)}^ , 

r 

where F goes around = 1. In particular, 

p(C) = ^^it*(C)(g^"^(°))}, 

They are given in terms of the eigenvalues of the transfer matrices and the Q op- 
erators, as well as other characteristics in the Matsubara direction. Their explicit 
formulas will be given in 02.61) and (17.21) . In Appendix C we shall explain that uj{(,C,) 
is a quantum deformation of the canonical normalised second kind differential on a 
hyperelliptic Riemann surface. 

From the equations (11.91) . (ll.lOp . (11.111) one immediately derives 

(1.12) Z^{t*iO ■ • •t*(4°)b*(Ci+) ■ • ■b*(C+)c*(Cr) ■ ■■c*iCi){q"''^'^)} 

= n2p(C°)xdet HCC-)),^.^,...,. 
p=i 

Taking the Taylor coefficients in [Q)'^ — 1 in both sides, one obtains the value of Z'^ 
on an arbitrary element of the basis (II. 8p . This is the main result of the paper. 

The text is organised as follows. 

In Section 2 we give the precise definition of the linear functional Z^ on the space 
Wa,o- We explain that on any particular X G V^a,o this functional reduces to a finite 
expression. 

In Section 3 we prove (II. 9p . A significant part of this section is devoted to the 
reduction of Z'^|t*((^)(X)} to finite intervals. This is a point which is used in Section 
6. 

In Section 4 we explain some simple facts about transfer matrices and Q operators 
in the Matsubara direction. It should be considered as preparation for the following 
sections. 

In Section 5 we introduce g-deformed Abelian integrals which are constructed via 
eigenvalues of Q operators in the Matsubara direction. We introduce g-deformed 
exact forms and present the g-deformed Riemann bilinear relations. 

In Section 6 we consider Z''|b*((^)(X)}. We formulate two lemmas which are 
proved in Appendix A and Appendix B. Informally, these lemmas say that Z'^|b*((^)(X)} 
is a g-deformation of a normalised second kind Abelian differential in (, which has a 
prescribed singularity specified by the quasi-local operator X. In the classical limit. 
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such a differential can be expressed using the canonical normalised second kind dif- 
ferential. Formula fll.lOp is an analogue in the quantum case, the function uj{(,C,) 
playing the role of the canonical differential. 

In Section 7 we define uj{(,^). Using the results in Section 5, we prove that it 
satisfies all the necessary requirements. 

Finally, in Section 8 we prove the main Theorem which states that (11.101) . (11. lip 
hold. 

As mentioned above. Appendices A, B are devoted to the proof of the technical 
Lemmas in Section 6. In Appendix C we consider the classical limit of the q- 
deformed Abelian integrals and differentials. Then we explain that the classical limit 
of uj{(,^) is indeed related to the canonical normalised second kind differential. Some 
general information about differentials on Riemann surfaces is provided. Readers 
who are not familiar with Riemann surfaces are recommended to read Section 5 and 
Appendix C together. In Appendix D we show equivalence of several non- degeneracy 
conditions accepted in the text. 

2. Definition of the linear functional 

Consider a two dimensional finite lattice composed of two one-dimensional chains: 
the space chain, and the imaginary time or the Matsubara chain. The space chain has 
21 sites which are labelled by the letters j = — / + 1, ■ ■ ■ , / . With every site the Pauli 
matrices a" are associated. The Matsubara chain has n sites labelled by boldface 
letters m = 1, ■ ■ ■ , n. With every site we associate a half- integral spin s^n and a 
parameter in other words a (2sjn -f l)-dimensional evaluation representation of 

the quantum group Ug{sl2). We assume that X]m=i ■^m is an integer. 
We define the monodromy matrix 

Tj-m(C) = ^j,n(C/^n)i^j,n~l(C/^n-l) " " " i^j,! (CM) • 

The L operator i^j,m(C/Tm) is obtained from the universal one 



^ [ Cq - q 2 [q-q ^)(Fq 2 
{q — q^^)Cq ~E (^q 2- — q^~ 
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by letting F, H act on the (2sin + l)-dimensional representation of Uq{sl2)- We 
shall consider a twisted transfer matrix 

Tm(C, i^) = Trj (rj-M(C, /«)), 

Tj,^{C,K)=Tj,M{Oq^''l 

and use the letter T{(, k) to denote its eigenvalues. 

Now we are ready to introduce the main object of our study. On the space 
consider the linear functional 

^ i l^'^ TrMTr[_,+i,,] (Th+1,/],M g2(«5[_,+,,,+a5[_,+,,o,)) 

Here and for later use, we set 

(2.2) T[k,m\,M = Tk,M ■ ■ ■ Tm,M , ^j,M = ^j,M(l) • 



6 



M. JIMBO, T. MIWA AND F. SMIRNOV 




q 



K a 



fig- 1 



In this picture the broken hnks represent the operator 0: the arrows on them are 
fixed. 

Suppose that the transfer matrix Tm(1, has a unique eigenvector \k) such that 
the corresponding eigenvalue T{1,k) has the maximal absolute value. Similarly let 
{K + a\ is be eigencovector of Tm(1, n + a) with the eigenvalue T(l, K + a) possessing 
the same property. Let us remark that for the XXX model the spectrum in spin 
zero sector is known to be simple even in the homogeneous case [13] even in the 
homogeneous case. Suppose also that 

(2.3) {K + a\K)^0. 

It is clear that in this situation (12. ip reduces to the linear functional 



(2.4) Z'^(g2"-5(o)o| = lim 



{k + a Tr[_i_ 




l+l,l],M Q ^ ' 


,+i,i]+aS[_,+i,oi)oj 


(/t + a Tr[_ 


+1,1] (^[ 


-«+i,/],M q ^ ' 


-i+i,i]+a5'[-i+i,o])^ 



This is the object which we shall calculate. For any given quasi-local operator we 
can proceed further. Indeed, if the support of q'^°'^^^^0 = q'^°'^^''~^^X[k^m] is contained 
in the interval [k, m] of the space chain, then 

ynj2aS{k-l)y \ „,,^fc-l (^ + «|Tr[fcH (%m],Mg'"^"'-"''X[fc,^]) \k) 

^1^ ^X[.,^])=P(1) rfl,«:)-^+i U + a\n) 

where 

(2.6) piO 



r(l,K)"^-'=+i {K + a\K) 
T(C,a + K,) 



The function (12.61) will play an important role for us; we shall see in the next Section 
that this is the same function as in (II. 9p . The last formula (12.50 shows, as it has 
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been said, that the hmit / ^ oo is superfluous. It is put in the formula (12. 4p just 
for the sake of treating all quasi-local operators simultaneously. 

It may look surprising that the thermodynamic limit in this approach is so sim- 
ple. Usually, it requires a complicated analysis of Bethe equations. Certainly, the 
complexity of the problem cannot disappear, and it is hidden in the limiting process 
n ^ oo to arrive at (11.31) . But the idea used in [7], and developed further in [8], is 
that one can proceed rather far before taking this limit. This is especially true in the 
present work. The complexities of the thermodynamic limit of Z'^{X) are confined 
to only two functions, for which one can take the limit n ^ oo rather easily. 

Let us emphasise one point which may be a source of confusion. We started with 
(12.11) . reduced it to (12. 4p and further to (12. 5p . The expression in the right hand 
side of ( 12. 4p is perfectly well defined for any pair of eigenvectors of Tm{Cj^) and 
^m(C; «^ + «) satisfying the condition (12. 3p . For the computation of (12.40 we shall use 
only quite general facts concerning Bethe vectors, so they are valid in general. Still, 
the subject of our study is (12.11) . and it reduces to (12.41) only for the eigenvectors 
corresponding to the maximal eigenvalues. 

3. Computation of Z'^{t*(C)(X)}. 
According to (II. 7p we are actually interested only in the following block of t*((^): 

t*(C,«)=t*(C)lw.,o-.W.,o- 

Let us recall the definition of the operator t*{(, a) given in the paper [2]. We start 
with a finite interval and an operator X[k,m] ■ With this notation we imply that X^k^^n] 
acts as / outside [fc, m]. Define for Z > m 

t[fc,«](C,«)(^[fc,m]) = Tr„ {Ta,[k,l]{C,a){X[k,m])) , 

where 

Ta,[fc,«](C; tt)(-^[/c,m]) = ^a,[fc,«] (C)Q'"'^"-^[fc,m]^a,[fc,/] (C)~\ 
Ta,[k,l]iC) = Ra,liC) ■ ■ ■ Ra,kiO ; 

Ra,j{C) is the standard 4 by 4 i?-matrix (see e.g. (2.4), [2]). Define further 

Rl^iC') = C'"M,,(C)P,,C"'^ = 1 + (C' - l)r,,(C') , 

where IPjj(-) = Pij{-)Pij and Pij is the permutation operator. Since IRi^j(l) = 1, 
Tiji^C^) is regular at = 1. Then [2] 

/-I 

tf,,,] (C, «) (X^kM) = 2 E - (C') ■ ■ • rmwi iO^y iC) (F[fc,^+i] ) 

j=m 

+ (C' - l)'-'"Tr„ {r„XC')rM-i(C') " ■ ■ r„wi(C')K^(C')(>^[fc,™+i])} , 

where y[fc,m+i] = (l°"^'^'''{X[k,m]), t is the shift by one site of the chain to the right, 
and 

M^(C')(F[,,^+i]) = C+i,„(C') ■ ■ ■ir+l,fe(C')(l^[fc,n.+l]). 
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Hence the limit / ^ oo is well defined as a power series in ("^ — 1: 

oo 
j=m 

We repeated these definitions because we want to make clear the following point. 
Take a 2 by 2 matrix K such that Tr(i^') 7^ and consider the following object: 

2 

Then it is easy to conclude from the above definition that 

(3.1) tf,_,](C,«,ir)(X[fcH) =tffc,,](C,«)(X[,,™]) mod (C'-l)'-™. 
Lemma 3.1. We have 

(3.2) Z"|t*(C)(g'"^(°)0)} = 2p(C)Z"|g'"^(°)o} . 

Proof. Without loss of generality, let = Xji ^] be localised on the interval [l,m]. 

(/c + a|Tr[i,;],jT[i,i],Mg'"''i^-"T,,[i,;](C,a)(X[i,™])]|/s:) 

= lim -7- — 

1^00 T(l, k)'(k + q;|k) 

From the considerations above we obtain 

(/c + a|Tr[i,i],jT[i,;],Mg'"''i^-'iT,,[i,i](C)(X[i,„])) |k) 



^ {k + a|Tr[i,;], Jr[i,z],Mg''^^'^-"r,,M(C)g''"'Ta,[i,^](C, «)(^[ih)) \^) 



mod (C^ — 1 



The idea here is exactly as in (13.11) . The monodromy matrix Ta;M.{.Cj<i'^"°^ plays 
the role of Ka- The fact that it carries the additional structure as operator in the 
Matsubara space is not important. What is important is that the state \k) is an 
eigenstate of Tra(Ta_M(C)Q''"^") with eigenvalue T{C,k). Now we can proceed using 
the Yang-Baxter equation and the cyclicity of trace: 

7f^{'^ + «|Tr[i,/],a(r[i,,],M g''^^^'iT„,M(C)g"'''Ta,[i,^](C, «)(^[ih)) I'^) 
= J^^^^ + «|Tr[i,,],,(T,,[i,,](C) (T,,M(C)g^'^+"^'^'T[i,,],M |k) 



nc/^) 

= 2p(C)(ft: + a|Tr[i,^] {Tyi,i]M \^) ' 

which proves the assertion. □ 



GRASSMANN STRUCTURE IN XXZ MODEL 



9 



Some comments on (13.21) have to be made. It has been said that r = t*/2 is the 
shift by one site of the chain to the right. According to [2j the rest of t* is constructed 
from the adjoint action of local integrals of motion. Then, looking at (12. ip one may 
wonder where p{() comes from. Naively, it should not be in the right hand side 
because r and adjoints of the integrals of motion commute with Ttm(Tq^m) and 
hence they should not contribute to (12. ip due to the cyclicity of trace. However, 
this is not correct in the presence of disorder field Let us explain this point 

in the simplest case r = t^/2. Consider the definition (12. ip . For finite / in (12. ip . we 
define the cyclic shift by one site t'p™"'^^'^ ^ which acts in particular as 

^periodic /^2aS[_i+i,0l A = g2o5[_;+2,l] 



On the other hand, it is easy to see from the definition that our operator r acts as 

This difference accounts for the appearance of p(l) in the functional (12. ip . A similar 
thing happens with the adjoint action of the local integral of motion Ip(-) = [Ip, ■]. 
The operator ]lP<5"°dic fgg[g ^^^q inhomogeneities of q'^°'^i-w,o] ; between sites and 
1 and between sites — / + 1 and /, while the operators entering the definition of t*{() 
feel only the first one. This is the reason why p{() appears. There are two cases 
when p{() = 1. The first one is trivial: a = 0. The second one is the case of VEVs 
(ll.2p which was considered in pQ, [2]. 

Before proceeding to b* and c* we have to give some explanation about g-deformed 
Abelian integrals. 

4. Spectral properties in Matsubara direction. 
Consider the transfer matrix 

TM(C,A) = Tr4r,,M(C)g''^"). 
Let us introduce the Q operator 

Qiti(C,A)=C'-^Tr^(TA,M(C,A)), 
where S is the total spin operator acting on the Matsubara chain, and 
Ta,m(C, A) = LAAC/rn) ■ ■ ■ LA,i(C/ri)g'"^^ . 

Here the L operators are associated with the g-oscillator algebra with generators 
ra, aL*A, Da- For the notation and conventions, see [2j. If Sm = 1/2, then 

'l_^2g2D^+2 

To obtain LA,m{() for other spins, one applies the standard fusion procedure. The 
Q operator (5m(C) A) is defined by 

Qm(C,A) = JQ+(C,-A) J, 

where J is the operator of spin reversal. 



(4.1) LaAO 
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These Q operators satisfy the Baxter equation: 

(4.2) Tm(C, A)Q± (C, A) = diOQUCq, A) + aic)QUCq~\ A) . 

These equations hold for the eigenvalues because Tm(C) A) commute with Q^{^, A). 

The functions a{(), d{() are defined by spins and inhomogeneities present in the 
Matsubara direction. 

n 

(4.3) a(C) = n ««^(C/rm), a.(C) = CV^' - 1 , 

m=l 
n 

d{c) = n ^«^(c/^m), 4(c) = eq-''^' - 1 ■ 

m=l 

Let us cite one formula from [5]: 

(4.4) Q+ (C, A)Qm(C9, A) - Q^iC, A)Q+ (Cg, A) = ^,.s 1 ^-x+s ^iO, 



where 

n 2s 

w{() = n ^s^i(/r^), ^sio = n(i - cv'-'-^') ■ 

m=l k=l 

Suppose that Tm(C) A) has a unique eigenvector |A) with eigenvalue T{(, A) such 
that T(l, A) has maximal absolute value. We denote by Q^{^, A) the eigenvalues of 

n 

Qm(^' A) on I A). If the eigenvector |A) has spin (i— ^ s^, it follows from the form 

m=l 

of the L operator (jH]) that C~^^^Q^(C, A) is a polynomial in ("^ of degree d, while 

n 

C^~^(5~(C) A) is of degree 2 J2 — d. Due to the quantum Wronskian relation 

m=l 

(14.41) . their leading and the lowest coefficients are both nonzero. 
Let us discuss the symmetry under negating A. We have 

(4.5) Tm(C,-a) = jrM(C,A) J, 

which implies that the spectra of Tm(C) A) and Tm{(, — A) coincide, and, in partic- 
ular, 

(4.6) T(C,A)=T(C,-A). 
Furthermore, the equation 

Qm(C,A) = JQ+(C,-A) J 

implies that 

(4.7) Q-(C,A)=Q+(C,-A). 

Due to (14.51) the vectors |A) and | — A) have opposite spins. 
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5. Deformed Abelian integrals 

Working with quantum integrable models, one should not neglect the important 
piece of intuition provided by the method of Separation of Variables (SoV) discovered 
by Sklyanin [13]. It has been explained in [TU] that the matrix elements of observables 
in the SoV method are expressed in terms of deformed Abelian integrals. In the 
case under consideration, which is related to the algebra Ug{sl2), these integrals are 
deformations of hyperelliptic ones. Let us give their definition. 

Introduce the function ip{() which satisfies the equation 

(5.1) aiCqMCq) = d{CMO ■ 

This function is elementary, 

n 2s ^ 

^iC) = n ^sAC/Tm) , VsiO = n ^2^-2s+2fc+l _ ^ " 
m=l fc=0 ^ ^ 

In addition to the contour F which encircles = 1, we consider n + 1 contours in the 
plane: Tq which goes around 0, and Fm which encircles the poles = r^g^'*'""^'^"^ 

(fc = 0, ■ ■ ■ ,2Sm) of (fs^iC/Tm)- 

In the following, we use the g-difference operator 



-1^ 



Ac/(C) = fiCq) - fiCq 

It acts on the class of functions of the form f{Q = (^^P{('^), P being a polynomial 
in and ^ 1 for all integers n. Within this class the g-primitive A^^/(C) is 

defined uniquely. 

There are two kinds of deformed Abelian integrals. 



(5.2) 



where C'^"/^(C) is a polynomial in (^^, in order that the integrand is single-valued. 

We start our study of deformed Abelian integrals with the following technical 
lemma. 

Lemma 5.1. Let (^"'f^iC) be a polynomial in ("^ . Then, for m = 0, ■ ■ ■ , n, the 

following identities hold. 

(5.3) 

I {Tie, «:)A^-V^(Cg) - Tie, K + «)A^-V±(C)}q^(C, k + a)Q^iC, nUQ^ 
= / fHOaiOQ^iCn + a)Q^iCq-\nMO^, 
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(5.4) 

= I /^(CKC)Q^(C,«: + «)Q^(Cg,'^MC)^. 

Proof. This can be verified directly by applying the Baxter equation to T{(, k)<5^(C) i^), 
T(C, K + a)Q^{(, K + a) and moving contours of integration. □ 

It is well known that, on a compact Riemann surface of genus g, the space of the 
first and the second kind differentials (meromorphic differentials without residues) 
has a finite dimension 2g when considered modulo exact forms. In [10] it was 
explained what are the g-deformed exact forms, for which the deformed Abelian 
integrals vanish. Since the proof was omitted in that paper we include it here. 

Lemma 5.2. Define a q-deformed exact form to be an expression 
(5.5) E{fHO) 

= T(C, k)A-i (/±(C)T(C, k)) + T(C, K + a)A-' (/±(C)T(C, + a)) 

- T(C, k)A^-i (/^(Cg)T(Cg, «: + «))- T(C, ^ + a)A^' {f^{Cq~')T{Cq~\ «:)) 

+ a{Cq)d{Of^{Cq) - d{Cq~')a{Of^{Cq-') , 

where C^"/^(C) « polynomial in C,"^ . Then we have 




Proof. Let us divide the integral into two pieces: 

(5.6) j E{f^iC))Q^{C, n + a)g±(C, ^MC)^ = h + h. 

where Ji first four terms from (15.51) and I2 contains remaining two. Apply (15.31) to 
the first and the fourth terms in /i, and to the second and the third terms as well. 
Using the Baxter equation and moving contours using (15. ip . one obtains 

h = j fHCq-')nCq-\'^)Q^{C,'^ + a)Q^iCq-\^)a{CMC)^ 
- I fHOT{C,K + a)Q^{C,K + a)Q^{Cq-\K)a{CMO^. 
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Now apply the Baxter equation: 

- / /^(C) {a{C)Q^{Cq-\^ + a) + d{OQ^{Cq,n + a))Q^{Cq~\K)a{CMO^ 

Moving contours we find 
h = J {a{Cq)d{Of^{Cq) - d{Cq-')a{OfHCq'')}Q^{C, n + a)Q±(C, ^MO^ , 

i.e. Ji = -l2. □ 

A beautiful feature of deformed Abelian integrals is that they allow for a defor- 
mation of the Riemann bilinear relations as well. In [10] the latter are given in the 
full-fledged form. For our present purposes, it is sufficient to use a part of them 
given by the following Lemma. 

Lemma 5.3. Consider the following function in two variables 
where 

r+(C, = r+(C, «), C) = r+(e, CI - -«), 

and 

(5.7) r+(C, a) = r(C, k)A^' m/^, «)(T(C, - m, ^ 
+ r(C, + a)A^' (^(C/e, a){T{C, + «) - r(e, + a))) 

- T(C, /t) A^-i mC/^, «)(T(Cg, «: + «)- T(e, + «))) 

- r(C, + a)A^' (^(g-'C/e, «)(r(Cg-\ - 

+ (a(Cg) - a(0)rf(C)^(gC/e, a) - KCg"') - d{0)a{Cmq-\/^, a) . 

Then 

(5.8) J yr(C,OQ"(C,«: + a)Q+(C,«:)Q+(e,/^ + a)Q"(e,«:)^(C)^(O^^ = 0. 
r- r- 

Proof. The proof is similar to that of the previous lemma. We apply Lemma 15.11 
invoke the Baxter equation and move the contours. When the Baxter equation is 
applied to expressions like ct)(^(C) '^) ~ k)) separately with respect to C 

and ^, a singularity may appear from ip{C/C,, «). In general, by moving the contours 
such a singularity produces intersection numbers as in the genuine Riemann bilinear 
relations (see pjj). In the present case this does not happen because the contours 
do not have nontrivial intersections. □ 
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Clearly C,"'r'^{C,^) is a polynomial in and ( C) is a polynomial in C^, 

both of degree n. This allows us to define the polynomials by 

n n 

m=0 m=0 

Introduce the (n + 1) x (n + 1) matrices 

(5.9) ^5 = I C^"+^jg^(C, n + a)Q^{C nMC)^ , 

(5.10) = I C^>f (C')g^(C, + «)Q^(C, «:)¥'(C)^ • 
Then (15 .Sp reads as 

(5.11) s+(yi-)* = yi+(3-)*. 

We explain in Appendix C that, in the classical limit g — > 1 (and for a = 0), A^, 
"B^ reduce to the matrices of a-periods of differentials of the first and the second 
kind, respectively. The relation ( 15. lip becomes one quarter of the classical Riemann 
bilinear relations which state that the full matrix of a- and 6-periods is an element 
of the symplectic group. 

Before closing this section let us make a comment. Suppose C^"/(C) is a rational 
function. We assume that the poles of this function do not overlap those of V'(C) 
and = 0. In this case, the g-primitive A^^/(^) is not uniquely defined, and in 
general develops infinitely many poles {n ^ Z, w are the poles of C^°'f{C) )■ 
Nevertheless Lemma 15.11 remains true. Actually it tells that the deformed Abelian 
integrals in the left hand side of (15. 3p . (15.41) do not depend on a particular choice of 
the g-primitive. For the same reason, deformed Abelian integrals of the g-exact form 
in Lemma 15.21 have unambiguous meaning. Later on we shall deal with examples of 
such g-primitives of the form A7^(?/'(C/^, a)P(C^)) or A7^(^(^/C, a)P(C^)). 



6. Properties of Z'={b*(C)(X)} and Z''{c*{()X)} . 

Our strategy is to compute Z''|b*(C) (X) } and Z''|c*(C)X) } inductively, reduc- 
ing them to similar quantities involving the annihilation operators h{(), c{(). It has 
been said in Introduction that Z''|b*(C) (X) } is non-trivial only when X G W^+i-i, 
and Z'*|c*((^) (X) } is non-trivial only when X E 'Wa-ii- We denote these blocks 
by 

b*(C,«) = b*(C)lw„+,._i^w..o ' c(C,«) = c(C)lw.+i,_i-.W.,o ' 
c*(C,a) = c*(C)lw._i,i^w„,o ' b(C,a) = b(C)lw„_i,i^w„,o ■ 
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Hence the non-trivial part of our main equations (11.101) . (11. lip take the form 

(6.1) Z^{h*{C. a){X)} = ^J uiC, OZ^{c{^, a) (X) } ^, X G W,+i,_i , 

r 

(6.2) Z^{c*{Ca){X)}=-^Jco{C,OZ^{H^,cy){X)}^, X G W, 



Our task is to estabhsh the existence of uj{(,^) and to determine it exphcitly. In view 
of the spin reversal symmetry which relates (b*, c) with (c*, b), we shall concentrate 
on the first pair. 

Apart from an overall power of (, h*{(,a) is defined a priori as a formal power 
series in — 1. Nevertheless when acting on each operator it reduces to a rational 
function, due to the same mechanism as explained for t*. Namely, 

(6.3) b*(C,a)(g2("+^)^(°)X[i,„]) = g^^^^ lim Tr,{T,,[„+i,,](C)g.,[iH(C, «)(^[ih)}- 

/ — >oc 

We recall the definition of the operator gcji,m](C) ck) in Appendix Rl The formula 
(16.31) together with the requirement of translational invariance can be considered as a 
definition of b*(^, a), but self-consistency of this definition requires that gc,[i,m](C) Q^) 
satisfies certain reduction relations which were proved in [2]. 
Using (16.31) we find by the same method as in Lemma 13.11 

(6.4) r(C, k)Z^ (b*(C, «)(g^("+^)^WX[iH)) 

~ T{1,k)'^{k + a\K) 

Due to this equation the left hand side happens to be up to the overall multiplier 
a rational function of ("^ with poles only at = q^"^. Its singular part is given 
as follows. 

Lemma 6.1. Set 

(6.5) a;,,„,(C,0 = -Ac^(C/e,a) 
Then we have 

(6.6) T(C,/s:)Z'^{(b*(C,a) - ^ j 0Js^nAC.^)^{^.o^)){X)]'^ = CP^iC), 

r 

where X G V^a+i-i, T encircles = 1, and -Pn(C^) « polynomial in ('^ of degree 
at most n. 

Lemma 16.11 is proved in Appendix |X1 

In order to characterise the quantity in the left hand side of (16.61) , we need to have 
a control over the unknown polynomial Pn(C^)- This is the point where deformed 
Abelian integrals come into play. Introduce the notation 

D^F{C) = F{qC) + F{q-\) - 2p(C)F(C). 
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Lemma 6.2. For m = 0, ■ ■ ■ , n, the following relations hold. 

(6.7) I T(C, n)Z-{ (b*(C, «) + ^ / ^ (DcD^A-V(C/e, «)) c(e, «)) (X) 



xg-(c,«: + a)g+(c,/t)^(C)^ = o, 



C 



As explained at the end of Section [5l one can apply Lemma 15.11 to f~^{C) = 
D^ipiC/^, a). Then, the integral over ("^ in the second term can be rewritten as 

(6.8) J T(C, k)D^D^A^'^{C/^, «)q-(c, k + a)g+(C, kMC)^ 

= j D^^{C/^,a)Q'{C,K + a) {a{C)Q+{q~\, ^) - d{OQ+{qC, k)) ^{Q^. 

Hence it does not actually depend on a particular choice of A'^'^ipiC/^y ct)- 
Proof of Lemma [6.21 is long and technical. We defer it to Appendix [Bl 
Comparing (16.11) with (16.61) . (16.71) . we infer that the function uj{(,^) = uj{(, a) 

satisfy the conditions 

1. Singular part 

(6.9) C^TiC, K.) (t^(C, - ^^singiC, 0) is a polynomial in C of degree n . 

2. Normalisation 
(6.10) 

j T(c, k) {uic + :Dc^cA^" V(C/e, «)) Q~(C, + «)Q+(C, f^MC)^ = o , 

(6.11) (m = 0,---,n). 

Furthermore, the equation (16.21) requires an additional property of Lj{(,C,\K,a) (see 
Section 8). 

3. Symmetry 

(6.12) o^(e,ci-'t,-a)=^(c,ei/^,«). 

7. Definition of uj{C,^) and its symmetry 

We shall first give the definition of the function uj{(,C,), and then prove that it 
satisfies all the necessary properties. 

In Section [5l we defined the matrices and B"*". In Appendix [D] we show that 
the condition 

(7.1) detyi+^0. 
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is equivalent to the non- degeneracy condition fl2.3p accepted previously. The classical 
analogue of fl7.ip states that "there are no holomorphic differentials such that all 
the a-periods vanish" . 

Assuming (17. II) . consider the function 

(7.2) u{C, a) = -^^-i-^^;+(C)*(yi+)-iS+t;-(0 + uj^yUC, ^l^^, «) , 

where are vectors with components = "B are given by (15.91) . 

fl510|) . and 

UsyUC ^\^, «) = (4a(0rf(C) - T(C, n)Ti^, n)) ^(gC/^, «) 

- (4a(C)rf(0 - T(C, «:)T(e, /«)) ^(g-^C/^, «) 

-2V^(C/e, «) (t(C, «:)r(e, /€ + «)- r(e, «:)r(c, «: + «))}. 

The function UgymiCy^l^yCi) satisfies the relation 

(7.3) iOsymiC, ^|/«, a) = UJsymi^, CI " 1^, 

due to (14. 6 p and the equality ^'(C^"'^) = ~V'(C) Q;)- 

The function C~"i^(C) 'Cl'^; ci;) is a rational function of C^. It is clear by the con- 
struction that the property (16.90 is satisfied. 

The remaining properties (I6.10p . (16.120 of uj{(,^\n,a) are more comphcated, and 
we formulate them as lemmas. 

Lemma 7.1. The function u!{(,C,\n, a) defined by (17. 2p satisfies the normalisation 
condition (16.100 . 

Proof. By using definitions (15.90 and (I5.1U0 we have, 
The definition (15.70 can be rewritten as 

r+(c,o = ^(^(c/e,«)) 

- iT(e, K)TiC, /t) {uJsymiC, ^l^^, «) + ^Df^^ A^" V (C/^, «) } • 

Therefore, T{(, k) {uj{C,,^) + Dc_D^^^^il){C,/ a)) is a g-deformed exact form in C,. 

□ 

Lemma 7.2. The function u!{C,^\k, a) defined by (17.20 satisfies the symmetry con- 
dition (167121) . 
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Proof. In Section [5] we had the relation (15.111) . 
(7.4) 'B+{A-y = A+{'B'y . 

In Appendix |D] we show that det{A~) ^ follows from the condition (12. 3p . Hence 
both can be inverted. So, invert them and multiply the result by v^[C^^ from 
the left and f~(0 from the right: 

t;+(C)* [A-^Y^ S+t;"(0 = v-{iy [A-y^ '^'v^Q . 

What remains to do is to add oOsymiC,-, ^l^t, «) to both sides, to use (17. 3p and to recall 
the identities (USD and (gZH- □ 

8. Main theorem 
Now we are able to prove our main theorem. 
Theorem 8.1. Under the generality condition (12. 3p we have 

(8.1) Z-{h\0{X)) = ^/c^(C,0^1c(0(X)}^, 

r 

(8.2) Z-{c*(C)(X)} = /^(e, C)^ib(0(X)}^ . 

r 

Proof. Consider (18.11) . It has been said that it is sufficient to consider the blocks 
h*{(,a), c(^,q;). Due to the structure of singularities (16. 6p and (16. 9p we have: 

(8.3) TiC,n)Z^{{h*iC,a) - -L ^ ..((, Oc(e, «)) (X)}^ = CP^iC') , 

r 

where -Pn(C^) is a polynomial of degree n. Due to Lemma [6.21 and Lemma [7.11 we 
have 

J C PniOQ' iC K + a)Q^ iC, kMO^ = 0, m = 0,--- ,n, 

which implies Pn(C^) = due to (O). 

Now consider (18. 2p . According to [,2j, the operators c*, b are related to b*, c by 
the transformation 

0„(x(C, a)) = q-^N{a - 1) o J o x(C, -a) o J, 
where N{x) = — and J(X) = JXJ~^ is the spin reversal. Namely, 

c*(C, a) = -0a(b*(C, a)), b(C, a) = 0a(c(C, a)) . 
It is also easy to see that 

Z^{X} = Z-^{I{X)}. 
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Hence (18.11) implies 

Z'^{c*{C,a){X)} = ^ uiCCl - ^,-a)Z^{hiC,a)iX)}^ , 

r 

which is equivalent to (18.21) due to fl6.12p . □ 

Appendix A. Proof of Lemma 16.11 

In this appendix, we prove Lemma 16.11 We also prove some additional result 
used in Appendix B (Corollary IA.2I) . First let us comment on the equation (16.40 . 
This formula is used in the proof in order to reduce the action of the operator 
b*, when it is considered inside the functional Z'^, to that of an operator on the 
interval [1,'ni]. This is a great simplification, because without Z'^ the support of 
the coefficients in the expansion of h*{(,a) in — 1 becomes indefinitely large. 
In other words, inside Z'^ the series expansion of h*{(,a) can be summed up to 
a rational function. Therefore, the proof of Lemma 16.11 consists in computing the 
singular part of the rational function. This task is done indirectly by considering 
an inhomogeneous space chain. We introduce inhomogeneity parameters $, = 
so that the original multiple poles ("^ = g^^ in the homogeneous chain are split into 
simple poles = q^'^C.j for 1 < j < m. Define the functional 



'■3 

T[l^m],Mi^y ^) — ^1,m(^1, k) ■ ■ ■ Tm,M(^m, ■ 

Using this functional the equation (16. 4p takes the form 
(A.2) Z''{b*(C,a)(g2(°+i)^WA[i,^])} 

= ^fi,m+l] {2gm+l,[l,m](^'"+l, «)(-^[l,m])} la = -=6n=l,Cm+l=C • 

First recall from [2j the definition of the operator k[i^m]((^, a) and its basic relations 
with C[i,^](C,a), C[i,^](C,a), i[i^m]{C,a)- 

(A.3) k[i,„](C,a)(X[i_^]) 

= Tr,,^ {(T+T|,,^},[i,^](C,«)C"-''^-K9"''''^-'^[i,m])} , 
(A.4) k[i,„](C,a)(A[i_^]) - A^f[i,„](C,a)(A[i_^]) 

l,m] l,m]) ~r l'[l,m] [l,m] ) • 

in the complex plane the operators C[i^m]{(, «), C[i,m](C) ct); f[i,m](C) ct) have singular- 
ities at ('^ = ^1 only. The operator gc,[i,m]{CjC() is given by 



(A.5) 2g,,[i,„](C,a)(A[i 

,m]) 

-2Tc,[l,m](C, «)f[l,m](C, «)(-^[1h) + 2Uc,[l,m](C, ")(-^[1h)' 
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where 



For the proof of Lemma [^TT] we compare the singularities of gc,[i,m](C) Q^)(-^[i,m]) 
inside the functional Z^^^^^^ with those of C[i a;)(X[i ,„]) inside ^[im]- 

It is known that gcji,m](C; ct)(-^[i,m]) is regular at = Now we compare 
res^2=^±2^2^gcji,m](C,a)(X[i,^]) with res^2=^2^C[i,m](C, a)(X[i,m]). Set 

Lemma A.l. The operator [a^, f/[i^m]]+ has its support in [1, m — 1]: 

(A. 6) [a^, U[i^rn]]+ = X[i^rn-l]Im, X[i^rn-i] = Tim {c^mU[l,m]) ■ 

We have 

(A.7) res^2=52^C[i,^](C,a)(X[i,„])— = f/[i,^]]+ , 

(A. 8) res^2=g2.^2^ (gcji.mKC. tt)(^[iH) + 1'c,[ih(C5 a)f[i,m](C5 a)(^[i,m])) 

Proof. Property flA.6l) appears in |3] as Lemma 2.6 (see also [2J, Appendix D). For- 
mula flA.7p is proved in [2J, Lemma 2.2. The calculation for flA.81) is similar, but we 
omit the details. □ 



Corollary A. 2. We have the relations between C[im]((^,a) and cyi^rn]{C: 



a] 



dC 

(A.9) res^2=^2(c[i,„](C,a) + t*i^„](C,a)c[i,„](C,a))(A[i,„])— = 0. 
Proof. To see flA.91) . it suffices to write 

reSf2=52^C[i^^](C,a)(X[i^m]) = Tr^ {cr~^'^a,ra^a\l,ra-l\{.im,Oi)U[l,m\) 
'^'^ ai^ a,ni^ a,[l,ni~ l](^m,tt)(^(0, l)mt^[l,m] (^q 



and use (lA.6p . (1A.7P and i?-matrix symmetry. □ 
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Using Lemma [A. II we obtain 



(A. 10) res^2^^^=,2g2^Zpi „+i] {2g„+iji,^](^^+i, 

= res^2=^2^Zfi„,] {c[i,^](C,a)(X[i,^])} — 

~ ^ ^6S«i+i=g2e^ri,»n+l] {f^[lH(^m^m+l " ^m^m+l)} ' 

Here the third term of flA.SP does not contribute, because the only singularities 
of i[i^m]{CyC() are the simple poles at = and the following inhomogeneous 
analogue of Theorem 13.11 holds: 

(A. 11) ^[i-m+i] {^rn+l,[l,m]{^m+l, Oi)X[i^rn]} = 2p(^m+l)^[Xm] • 

Note that 

r+(T+ — IT+T+ = (r^a^ — (T+t+ 

m m+l m m+1 \ m m+1 m m+1/ 

where P~ is the projector on the singlet, 

Using the cyclicity of trace and the quantum determinant relation 
we find 

m I 

^ Tr[l,m+l](/t + a|r[i,^_l],M($,/t)^[lH(^mQ^m+l - ^m^m+l)l^) 

{k + a\K,) 

-res^2=^2^Zfi ,„] |c[i,„](C,a)(X[i,„])) — . 



In the last line we used ( 1A.6p . (]A.7p . Computation of the residue at = q^'^i^ 



done similarly, using 

m+1 m m+1 m,m+l\ m m+1 ^m m+ll' 

The residues at C^+i = ^'^^ readily found from i?-matrix symmetry. 

Lemma A. 3. 

■sm+1 



(A.12) res^2^^^=g±2^2Zpi^^^_^i] {2g„+iji,^](^^+i, a)(X[i,^])} 



Sm+l 



res^2=g±2^2 u;(C,^j) res^2=g2Z['i {c[i,„](C, a)(X[i,^])} — . 
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Proof. This follows from the preceding calculations and 

reS(^2=^2^2 uj{C,Q = 1 - 



T{^j,K)T{^jq, k) ' 



□ 



Let us return to the homogeneous case = ■ ■ ■ = 6m = 1- The operators c(^, a), 
c{(, a) acting from Wq,+i _i to Wq, o are defined by 

(A.13) c(C,«)(g2("+i)^WX[iH) = g'"^(°)c[iH(C,«)(X[iH) , 

and the requirement of translational invariance. This definition is equivalent to the 
one given in [2]. The reduction relation proven there ensures the self consistency 
of the present definition. The equation (16.61) follows by writing (]A.12p as a contour 
integral and specialising to = ■ ■ ■ = = 1- 

It remains to show that the polynomial Pn(C^) iii the remainder term of (16.61) has 
degree at most n. The only non-trivial case to consider is when the spin of X^i ^] 
equals —1. Then it follows from the fact [2] that 



,[i,^](C,«)(X[i,„]) = o(i), c 



2 



OO . 



This finishes the proof of Lemma 16.11 . 

Appendix B. Proof of Lemma 16.21 

The goal of this section is to prove Lemma 16721 The proof is done in several steps. 
Step 1. 

Recall the definition (lA.Sp . Fix a solution foji,m](C,«) of the equation 
(B.l) Acfo,[i,H(C,a)(X[i 

,m]) k[l,m] 

which has poles only at ('^ = g^" (n G Z). Define further 
(B.2) b*(C,a)(g2("+^)^WX[i,„]) 

= lim g2"^(°)Tr,{T,,[„+i,i](C)go,c,[iH(C,«)(^[i,m])}, 

/ — >oo 

(B.3) go,c,[lH(C,")(^[l,m]) 

- Tc,[l,m](C; a)fo,[l,m](C; ")(-^[l,m]) + Uc,[l,m](C; a)(-^[l,m]) • 

Lemma B.l. Define 

= QmiC^ ^ + ")Tr[l,m],c(T'[i,m],M(l, k)Tc,m(C, (C,a)(X[i,„]))g+ (C,k). 
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The identity (16. 71) follows from 

(B.4) j /m(C)(X[i,^])¥'(C)^ = 0. 

Proof. Introduce the operator 

D^F{0 = F{qO + F{q-'0 - t*(C)F(C), 
which can be used to rewrite the definition of h*{(, a): 

+ hm Tre{Te,[„+i,,](C)ue,[iH(C, a)(^[iH)} • 

In similar formula for hQ{(, a) the only change is f to fo, u remains the same. 
Comparing this with the equations (1A.4I) . (1B.1I) and (1A.13I) we arrive at 

b*(C, a) - b*(C, a) = D^A^' (c(Cg, a) + ciCq'\ a) + c(C, «)) . 

Now consider the term containing c in (16.71) . Notice that due to Lemma [3.11 for 
any quasi- local operator X{() depending on ( 

Z^{D^{X{0)}=Z^{D^{X{0)}. 

So we replace in (16.71) D^, by D^^, D^. On the other hand from flA.9p it follows 
that we have an equality of formal power series in {(^ — 

(B.5) c(C,«)(g^("+^)"(°)X[iH) 

^ ^(C/e, «)t*(e, «)c(e, a)(g^("+^)^(°)X[,,„])^ . 



2Txi 

Using (IB.Sp we evaluate 



c(Cg, a) + c(Cg-\ a) + c(C, a) = ^ (^(C/^, «)) c(e, a)^ 



2 



In other words we obtain an equation similar to (16. 4p with b* replaced by the 
expression under in (16.71) . and g replaced by go, which is nothing but the matrix 
element of (lB.4p . 

□ 

Step 2. 

The next step is to reduce the identity to a difference equation for go on a finite 
interval. We will show that, for all m = 1, ■ ■ ■ , n, the identity (IB. 40 reduces to the 
same equation for a quantity in the space direction. So, we can forget the Matsubara 
direction. Introduce an operator 

[l,in] 

where 6 signifies the anti-involution 

e{x) = a'^x'a^ {x e End(r)). 
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Lemma B.2. Identity flB.4p follows from the equation 

(B.6) go,c,[lH(C,a)(^[l,m]) = -^c,[l,m](C) (gO,c,[l,m](^"^C,")(^[l,m])) " 

Proof. By symmetry it suffices to consider the case m = n. We prove the assertion 
assuming that Sn = 1/2. The general case is reduced to this case by the standard 
fusion procedure. 

From the defining relation 05.11] for ip{C), we have 

dC a{Tn) i^.dC 
resc2=g-2,2^(C)^ = ^^_^res^2=,2y^(C)^. 

So, the equation (IB. 41) is equivalent to: 

(B.7) d{T^q-^) /M(rn)(X[i,„]) + a(r„) Im{r^q~^){X^i^^^) = . 

Let us compute /M(Tn)(-^[i,jn])- We simplify notations introducing 

^) §0,c,[l,m 

First, move Tc^mi'Tn) to the left using the Yang-Baxter equation: 

= Tr[i,m],c(^rc,M(Tn, 't)7'[l,m],M(l, 't)7'c,[l,m] (Tii) ""^^[l,m] (Tii, a)^^!!,™] (Tji) j • 

Now, for Sn = 1/2, the L operator satisfies 

^c,n(l) = VPc,n, 

where we have set r] = q^/'^{q — q^^). Note that 

3 

where M' signifies the subinterval [1, n — 1]. Moreover, 

Tm{rn, K + a)Q^{Tn, k + a) = a{Tn)Q^{Tr,q~'^ , k + a) , 
because d{Tn) = 0. So we can evaluate /M(Tn)(-^[i,m]) as 

lM{Tn){X[i^rn]) = a(^n)QM (^ng"\ K + a) 

Now notice that 

^[l,m],M(l, 1^) = 7'[i,m]_n(Tn ^)^[l,m],M'(l) = /^(''"n)^n,[l,m] (^"n) ""^^[l.ml.M' (1 ) ) 

where /x(t) is a function whose explicit form is irrelevant for our calculation. 

Bring the permutation through T[i^m],M and put Tii^m],c to the right by cyclicity 
of trace: 

/M(l"n)(-^[l,m]) = Ai(^n)a(l"n)QM(^ng"\ K + a) 
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We compute /M(TnQ'^^)(-^[i,m]) similarly, using 
The result is 

X Tr[i,„],c(^T[i,m],M'(l5^)-Pc,n6'c(>ll,m](Tng~\a)6'c {Tc,[l,m]iTny^))^Q^{Tn, n) , 

where we applied the anti-automorphism 6c under Tr^ using Tr(^(x)) = Tr(x) and 
^c(2P-„) = Pc,n- Thus (El]) implies (ET]). □ 



Step 3. 

The third step is to reduce flB.6P to representation theory. 



Lemma B.3. Set 

yi = ^c^l - ^a(^t , 

Then equation ( ]B.6P ^5 equivalent to the identities 

(B.8) Trc,a,A (l/Tc,[l,m](C, a)Ta,[l,m](g~^C, «)TA,[l,m](g"^C, = 

{y = Z/1,1/2), 

where we have set = . 

Proof. Recall the definitions ( IB.ip . ( IB. 31) . From the definition of Ac^ji^m] we easily 
find that 

(B.9) Ac,[l,m](C)(XcX[i,m]) = Tc,li^rn]{C, a){X[i^ra])Oc{Xc) , 

(B.io) Acji,m](C)Tcji,m](g~^C,")(^[i,m]) = ^[iH • 

Write 

Applying (lB.9p . (IB.lOp we reduce (1B.6P to the form 

= (k[i,^](C,a)(X[i,^]) - Tc,[i,„](C,a)(k[i,^](g"^C,«)(^[i,m]))) 
+ (l[i,m](C,a)(^[i,m]) - Tc,[iH(C,a)(l[iH(^'^C,")(^[i,m]))) 

We rewrite this further, using the fact that ?/c = if and only if Tic{xcyc) = for any 
Xc e End(K;)- Nontrivial conditions arise from the choices Xc = or a~, giving 
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respectively 

(B.ll) k[iH(C,«)(X[iH) = Tr,,,,^((a+r- + {a^ - a^a+)a+) 

xT,,[i,„](C,a)T{,,A|,[i,^](g-iC,a)(g"'C)""''^""Kg~''''^""'^[iH))' 
(B.12) l[i,„](C,a)(X[i,„]) = Tr,,,,^((a+a- + {a^ - a^a+)r+) 



xT,ji,^](C,a)T{,,A},[iH(?''C,a)(g"'C)""^'^'-Kr''''^''"'^[iH)j- 
On the other hand, we have an identity (see pj, (2.20)) 

TA,[i,m](C,a)(^[i,m]) = Tr, (r+T|,,^},[i,^] (g-^C, «)(g"'C)°"''(^[iH)) ' 
which allows us to rewrite the left hand side of fIB.lip as 

k[l,m](C) Ol){X[l,m]) 

= TT,,a,A (a+r+TeJl,„KC,«)T{a,A}Jl,„^](g-'C,«)(g-'C)"-''^'™Kg"'^'^''"'^[lH)) • 

For l[i^m](C) '^)(-^[i,m]), we obtain an analogous expression, replacing o"^r+ by 
(a3 + a^a+)r+. 

After this rewriting we take the difference of the left and the right hand sides of 
(IB. lip . (]B.12p . and do further the gauge transformation 

Fa,A ■ T{a,A},[lH(C> «) " K,A = Ta,[l,m](C> a)TA,[l,m](C, ") (^a,A = 1 " ^A(^a) ■ 

The assertion of Lemma follows. □ 

To finish, it remains to prove (IB.Sp . Let "Jl be the universal R matrix of Uq{si2)- 
Denote by vr^ the evaluation module over V = C"^, and by ro^ that of g-oscillator 
representation W. For the notation and details, we refer to Appendix A of Let 
further vrji^] = vrf™". The identities of Lemma [B. 31 can be written as 

(B.13) Trc,a,A {y (ni^ Tig-ii; IS) rxjq^i^ nii^rn])^) =0 (y = ya)- 

In the tensor product 

Z = V^0 l/g-l^ ® Wg-ll^, 

there are two pairs which allow for non-trivial t/gp-submodules: 

Wo C Vg-l^ ® 

The submodule Vq — C (resp. Wq ~ VFg-2^) is spanned by f + (S> f- — f- ® f + (resp. 
{v. ® \n) + (g2" - 1)^;^ ® |n - l)}„>o). Set 

Then a direct calculation shows that 

yiZ d Zi, yiZi = 0, 

y2Z C Zi® Z2, 1/2^1 C Z2, 1/2^2 C Zi . 
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Since xZi C Zi {i = 1,2) holds for x G UqP, we have 

Trc,a,A (z/j(vr^ ® VTg-ic ® tZ7q-i^)(x)) = (x G f/^p) . 
The proof is now complete. 
Step 4. 

To complete the proof of Lemma 16.21 , we show matrix element of flB.4p between 
{k, + a\, \k) for m = 0. The integral consists of two parts, say Ji and J2, coming 
from the first 3 terms in (IB.Sp . or from Uc^[i^rn]{C , d) , respectively. We note that in 
view of (IB.ip and Lemma 15.11 the proper meaning of Ji is 



-^1 = ^ y + "|Tr[i,m],c(^T[i^„]^M(l,/t)^c,M(C,/«)kc,[iH(C5a)(-^[ 

To 

>-2 



xQ-{C, K + a) (a(C)Q+(g-^C, f^) - d{C)Q^{qC, n)) ^(C) 



The functions 

are all regular at ("^ = 0. On the other hand, for X[i^m] of spin —1, we have the 
estimate 

r"k[i,„](C,a)(X[iH) = OiC) iC - 0), 
as explained in p], section 2.5. The same argument there shows also 

C-"1[i,^](C,«)(X[ih) = 0(C') (C'-o). 

It follows that in both Ji and J2 the integrand is regular and the residue at = 
vanishes. This completes the proof of Lemma | 



Appendix C. Classical limit 

In this Appendix, we explain the classical limit of our construction and its relation 
to hyperelliptic Riemann surfaces. We shall not go into much details since similar 
considerations were done in [TU], [12] ■ We assume a = 0. At the moment we are not 
ready to discuss the classical limit in the case a ^ 0. We consider Bethe vectors of 
spin 0, so that C^'^Q^(C) ^) are polynomials in of the same degree 

n 

■S = ^ ■ 
m=l 

In the parametrisation q = e'^^'^, the classical limit amounts to z/ ^ 0. So z/ plays 
the role of Planck's constant. Let us see what happens to the solutions to the Baxter 
equation in this limit. First of all, in order to obtain a Riemann surface of a finite 
genus, we keep n finite. But for the classical limit we need to have large quantum 
numbers. This is achieved by considering large spins Sm- Actually, this was the 
main reason for us to consider arbitrary spins in the Matsubara direction. So, we 
require that usm, or equivalently g*"", tend to fixed non-zero values when z/ — 0. 
Similarly, we demand that uk, or q'^, stays finite in the limit. In this situation, a{(), 
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d{() and T{Q = T{(, k) tend to polynomials in (^^ of degree n, which we denote by 
the same letters as in the quantum case. 

In the classical limit, the poles of ^{Q concentrate on the portion of circles between 
the end points r^Q'"^*'" and r^q^'^"^. The s zeros of the polynomials C^'^Q^ICj'^) 
concentrate to n open curved segments Cm close to the above circular segments. 
This claim is difficult to prove, but we can justify them by analysing the Baxter 
equation in the classical limit. 

Consider the Baxter equation 

(C.l) d{QQ{Cq) + a(C)Q(Cg-') = T(C)g(C) ■ 

We look for its quasi-classical solution in the form 

(C.2) ^^^^ =^^^'^^^"^^2;?^ /logr/(0^}, 

1 

where 77(C) is a function independent of v and v) is a power series in u. First, 
dividing the Baxter equation by Q(C) and considering the leading order in Planck's 
constant we conclude that rj{Q must solve the equation 

(C.3) d{(uc) + a{c)r\c) = no . 

This is the equation of the spectral curve of the corresponding classical model. 
The function ri{Q has two branches, 

±^,. ^ no ± vncy - HOdjo 

^ 2d{0 
for future convenience we choose the branch of the square root such that a/ {q'^ — q'O"^ 

Consider the behaviour O' ^ 00. The polynomial T{0 is not arbitrary, it comes 
from the quasi-classical limit of a solution to the Baxter equation ( IC.ll) . Recall that 
for large C the function Q^(C) is O^C,^'^^"^") as discussed in Section 4. Also we have 
for (^^ — i> 00 and to the main order in Planck's constant 

a(C) = r-^g^^C'" + " " " , ^^(C) = r-^g-^^C'" + • • • , 
where r = H'^'m- So, the Baxter equation implies that 

T(C) = r-2(g'' + g-)C2" + ... . 

Hence when ^ 00 we have r/^ q±K+2s ^ which means that t]^ (resp. r/^) 
corresponds to quasi- classical limit of (resp. Q~). 

Throughout this paper we use as parameter C, while all the important quantities 
are actually functions of 0- This notational problem is due to historical reasons, 
and we are forced to tolerate it in the quantum case. However, in the classical case 
this notation becomes very unnatural making incomprehensible simple formulae for 
differentials on hyperelliptic Riemann surface. That is why in what follows we shall 
often use the parameter z = 0- For the same reason we denote the discriminant 
T(C)^ — 4a(C)rf(C), which actually depends on C^, by P(C^)- Recalling that a(C), 
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d{(), T{() are polynomials of and making the change of variables: z 
w = 2d{Qri{Q — T{() we bring the spectral curve fIC.Sp to canonical form: 



(C.4) = P{z) . 

In the g-deformed Abelian integrals the integration measure contains Q~{C, k)Q'^{(, k) 
The most direct way to compute this quantity uses the quantum Wronskian fl4.4p 

qK _ q-K 

^Av--v^)Q^iC,^)Q-iC,^)- 

This implies 

(C.5) Q+(C,/s:)g-(C,«:)^(C) 



.0 g---g-/p(C^' 

where we used the identity W{C)d{()ip{() = 1. 

The discriminant P{z) is a polynomial of degree 2n. Let us call its zeros Xi, ■ ■ ■ , X2n- 
The Riemann surface (1C.4P is presented as two copies of the z plane glued together 
along the cuts [a;2m-i, 2;2m]- According to the conjecture accepted previously, the 
branch points can be ordered in such a way that the cut [x2m-i, a;2m] is not far from 
the location of poles of v^(C) which are contained in the contour Fm- According to 
(IC.Sp . for any polynomial L(C^) we have in the classical limit 

(C.6) /l(CWC,«:)Q-(C,«:MC)-'^^' ^ f Hz) dz 



where Cm is a contour going in 2;-plane around [a;2m-i) a^2m] for 1 < j < n, or around 
for m = 0. The limit (1C.6P requires some remarks. The integral in the left hand 
side is taken over the contour T^. In the limit the integrand develops cuts which 
appear as a result of concentration of zeros of Q^{Ci ^)Q~(C) i^) a^iid poles of v^(C)- 
So, obviously, in the limiting process we have to deform the contour in order that it 
does not cross the cut. This is how the integral around Cm appears. 

The Riemann surface flC.4p has genus n — 1. The contours Cm, with m = 
1, ■ ■ ■ , n — 1 can be taken as a-cycles. Our Riemann surface have two points 0^ 
which lies on different sheets and project to 2; = 0. The contour cq goes around O"*". 
Similarly we have two points 00^ which project to z = 00. 

Define the differentials on the Riemann surface 

(j\{z) = — dz, j = 0, ■ ■ ■ , n . 

' ^/Piz) 

The differentials (Jj{z) where j = I,-- - ,n — 1, are holomorphic (the first kind) 
differentials while the differentials do and are the third kind differentials. The 
differential cxo has simple poles at z = 0^, it is dual to the contour cq. The differential 
CTji has simple poles ai z = 00^. 
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The holomorphic differentials can be normalised with respect to Ci, i = 1, ■ ■ ■ , n 

because 



det( / aj ^ 

i,j=l, ■■ ,n-l 



Ci 



This is the classical version of (17.11) . 

Consider the differentials whose only singularities are at 00=*=. Among those are 
exact forms 

(C.7) (^z^^/P{zj'^ dz, z^dz, k>0. 

Up to exact forms, holomorphic forms and the third kind differential (Jn{z) there are 
n — 1 linearly independent second kind differentials with singularities at cx)^: 



aj[z) = z- 



dz . 

J = I, - - ,n- 1 



2JP(z 



where [/(^)]+ means the polynomial part of f{z), which is a Laurent polynomial at 
z = 00. We shall use at some point the differential do which is an exact form. 

The most important identity in the theory of Riemann surfaces is the Riemann 
bilinear relations. Usually this identity is written in the form 

^ uji UJ2- uji uj2^ = 2ni U1OU2, 

m=l J K 

where lji 2 are the first or the second kind differentials, and 



UJi O UJ2 



^^res(co'i(i '^102)- 



In our case the a-cycles coincide with ci, ■ ■ ■ , Cn-i. The 6-cycle bm (m = 1, • " ' ; ^ ~ 1) 
crosses the cycle once on the first sheet of the surface, goes to the second sheet 
through the m-th cut, arrives to n-th cut by the second sheet, crosses this cut and 
returns by the first sheet to its beginning. 

An alternative way of writing the Riemann bilinear relations is the following. It 
is easy to see that a and a constitute a canonical basis 

(C.8) (Ji O CTj = 6ij, (Ji O (Tj = 0, CTi O CTj = . 

Now construct the antisymmetric form 

n-l 

(C.9) (T{x,y) = 5^(fTj(x)aj(y) - (Tj{y)aj{x)) . 

j=i 

Then 

(C.IO) j j aix,y) = 2mg^og^, 

gi 92 
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where in the right hand side we put the intersection number of cycles. From the 
exphcit formulae for and aj, one easily finds the 2-form a{x,y), 

This form is exact, so, apparently the integrals over all 2-cycles must vanish. How- 
ever, there is a singularity at x = y which produces the intersection number in the 
right hand side of (IC.lOp . All that is quite standard, so we do not go into much 
details. 

Consider a particular case of (IC.lOp . 



(C.12) / aix,y) = 0, i, j = 1, ■ ■ ■ , n - 1 . 




This is true because the a-cycles do not intersect. 

On the product of two copies of Riemann surface we have the canonical second 
kind differential p{x, y) with the following properties. 

• The differential p{x, y) is holomorphic everywhere except the diagonal, where 
it has a double pole with no residue 

(C.13) p(x, y) = r ^ + 0(1)^ dxdy . 



• The differential p{x, y) is normalised with respect to x, 
(C.14) Jp{x,y) = 0, m=l,---,n-l. 

Cm 

An important consequence of the Riemann bilinear relations is that this differential 
is automatically symmetric: 

(C.15) p{x,y) = p{y,x) . 

Let us explain this by giving an explicit construction of p{x,y). We start with an 
exact form in x, 

_d_ ( y/pR \ 
dx \^P^)[^j: - y) j 

which obviously has the required singularity at a; = but has also additional 
singularities at infinity. Because of flC.9l) and fIC.lip . these singularities are cancelled 
in the following expression: 

p(x, y) = -— 1^-^^=^— — y j dxdy + ^ o;{x)a;{y) + X;-ja^{x)a;{y) , 



dxdy . 



where the matrix Xi j must be defined from the normalisation condition 
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Now writing a similar formula for p{y,x), it becomes apparent that the symmetry 
flC.lSp is equivalent to the fact that X is a symmetric matrix. This fact follows 
from flC.12p . There is an obvious similarity between this argument and the proof of 
Lemma I7.2[ 

Suppose that we want to construct a normalised second kind differential with 
given singular part. To be more precise, we allow a singularity only at x = 1 with 
a given singular part 

N 
k=2 

So, we look for a differential which has the singular part Tsing{x) at x = 1 and is 
holomorphic elsewhere. We require that t{x) is normalized 

r(x) = . 

It is rather obvious that t{x) is given by 

(C.16) t{x) = (h a{x,y) d'^Tsing{y) , 



r 

where the contour F is as usual: 1 is inside it, and x outside. 

Let us return to the quasi-classical limit of the quantum formulae. First, notice 
that for a = the operator D becomes the second difference derivative because 
P(C) = 1: 

^^c(/(C)) = ^ (/(Cg) + fiCq-') - 2/(0) - (c^) /(O . 

Also A^^ goes to the primitive function 



2--V(/(C)) -„^C^ ) /(O 



Consider the /(O = -^(C^) ^-iid the corresponding g-deformed exact form (for a = 
there is no difference between /^(C)): 

^AC) = —E{m)Q-{OQHCMO% 

then 

^u(z) -j-^{L{z)^P{z))dz. 

Denote 

1 dC^ df"^ 

<yu{C\e) = —riCOQ-iOQHCMOQ-iOQHOviO^^, 

Then we have 

(C.17) (Ju{x,y) ^ (T{x,y) + ^{(7o{x)ao{y) - ao{y)&o{x)) , 



GRASSMANN STRUCTURE IN XXZ MODEL 



33 



the additional term is not important in fIC.lOp because ctq is an exact form. The 
hmit fIC.lTp explains the name g-deformed Riemann bilinear relations for fl5.8p . 
Consider now 

1 /7/-2 2 

p.ice) = —nomMcoQ-iOQHcMOQ-mHo^io^-^- 

We want to show that 

First, it is rather easy to find that in the singularity (16.50 two simple poles produce 
in the classical limit the double pole in (IC.lSp . Second, we have the normalisation 
conditions fl6.10p . They look different from the normalisation conditions (1C.14P 
because of presence of the term 

(C.18) — / T{c,^)D^D^A^'^{c/OQ-iOQ^{CMO%- 

However, this term for a = 0, z/ — is of order u"^, while PuiC'^,^'^) is of order 1. 
So the term ( ]C.18P does not count and from ( ]6.10p with m=l,---,n — Iwe get 



the normalisation conditions flC.14p . Conditions fl6.10p with m = 0, n show that the 
differential PuiC'^,^'^) in the limit u ^ does not have simple poles at = 0, cxo 
which were originally present. 

Thus we conclude that the function uj{(,C,) is related in the classical limit to the 
canonical normalized second kind differential. 

Notice a clear similarity between the formula (1C.16P and our main formula (16. ip . 

Appendix D. Equivalence of different non-degeneracy conditions. 

In this Appendix we show that the conditions det(yi^) 7^ are equivalent to the 
fact that the scalar product (12. 3p does not vanish. We use usual notations of the 
Quantum Inverse Scattering Method (QISM) 



//1(C) BiC) 
IaM[(.) - [^dO D{C) 

Consider the case when all the spaces in Matsubara direction are two-dimensional 
(spin 1/2). The basis of the two-dimensional space will be denoted by e±. Introduce 
two vectors in Matsubara space 

(D.l) 1+) = e+(g)---(g)e+, , |-) = e_ (g) ■ ■ ■ (g) e_ . 

The eigenvector |k) is written in QISM framework as 

(D.2) i'^)=n^(v)i-)' 

where (Aj~)^ are zeros of C'^Qm(C? /*) which is a polynomial of C^. It is well-known that 
this eigenvector does not vanish identically unless Ti = rjg for some j > i. The latter 
situation has to be forbidden from the very beginning because the tensor product 
on i-th and j-th spaces is reducible and contains one-dimensional sub-module. On 
the other hand there is no problem with the case Ti = rjg"^ which allows the fusion 
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procedure, and show that our considering only spin 1/2 representations is not a real 
restriction. 

Consider now the vector Y[ -B(Aj") |+), where (Aj*")^ are zeros of C~'^Qm k(C)- This 
vector also does not vanish identically, it is an eigenvector of Tj^{(,k) with the 
same eigenvalue as flD.2p . Hence, the assumed uniqueness of the eigenvector with 
the eigenvalue of maximal absolute value implies that this vector is proportional to 
\k) with some coefficient which depends on tj and k, the exact form of this coefficient 
is irrelevant here. 

Now consider the scalar product (12.31) . We do not care about the normalisation 
of the eigenvectors, so, in traditional QISM way it is written as 

('^+«i'^)=(-in^(/^7)nc'(V)i-)' 

where {^J'f)'^ are zeros 

of C^^QmIC^;)- Due to the previous remark we rewrite: 
(D.3) {K + a\K) = Const - {-\Y[B{fi7)YlB{Xf)\+) , 

where Const is a nonvanishing constant which was discussed above. So, we conclude 
that the scalar product in question is given essentially by the partition function with 
domain wall boundary conditions 

n 

M„(ei,--- ,en|ri,--- ,r„) = J]^ri(-| n^(^j) 1+) 

j=i 

with specification = {yUj^} U {Aj^}, notice that independently of spin of our 
eigenvectors the number of elements in the latter set is n. 

Being a polynomial of degree n — 1 in the function is completely charac- 
terised by the recurrence relation: 

(D.4) 

M„(ei,---,e n— 1 ; '^n I ''"l ) ■ ■ ■ ''"n— 1 > ''"n ) 

= (q' - l)^n ' n^^'n^^'^f - ^nKQ'rl - rj2)M„_i(ei, ■ ■ ■ ,en-i|ri, • ■■T^^^) . 
This recurrence was solved by Izergin who found a determinant formula for Mn_i 

m- 

On the other hand we have the determinant det(yi"^) of (n + 1) x (n + 1) matrix. 
This determinant depends on the Bethe roots only through the product Q~{C, k, + 
a)Q~^{(, k). Once again we consider the union = U {Aj"} and normalise 

this product as follows 

n 

j=l 

The determinant can be reduced in two steps: 
(D.5) det(yiJ)i,j=o,.,n = -27r2 J]e^det(yijj)ij=i,.,„, 

det(/ljj)ij=i,...,n = -27ridet(yiJ)ij=i,...,n--i , 
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where we used the obvious identities: 



Too 

Making the dependence on n and other parameters exphcit we introduce 

where we preferred the intermediate reduction from (ID.SP for its antisymmetry with 
respect to permutation of r's. In the case of two-dimensional representations in 
Matsubara direction the integrals in A^^ are easy: they are given by sum of two 
residues. Obviously, is a polynomial in of degree n. However the second 
relation from (ID.Sp shows that the actual degree is n — 1. 

Set = Tn and multiply the matrix from the right by the matrix I — t^E with 
Eij = (5j .j-i. Then it is easy to see that in the last row only n-th matrix element 
does not vanish. Using this, after some simple algebra one sees that satisfies the 
relation flD.4p . Hence we conclude that 

Due to the above reasoning it shows that {k + a\K,) is proportional to det(yi"'") with 
non- vanishing coefficient. Similarly, rewriting {k + a\n) as 

{k + a\K) = Const ■ {+\l[Cifif)l[CiXf)\-) , 

one proves that it is proportional to det(yi^) with non- vanishing coefficient. 
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